Neutron scattering and superconducting order parameter in YBa 2 Cu 3 07 
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We discuss the origin of the neutron scattering peak at 41 meV observed in YBa2Cua07 below T c . 
The peak may occur due to spin-flip electron excitations across the superconducting gap which are 
enhanced by the antiferromagnetic interaction between Cu spins. In this picture, the experiment is 
most naturally explained if the superconducting order parameter has s-wave symmetry and opposite 
signs in the bonding and antibonding electron bands formed within a CU2O4 bilayer. 

PACS numbers: 74.72.Bk, 61.12.Bt, 74.20.Mn 

Neutron scattering is an important tool in the study 
of high-T c superconductors. A sharp magnetic neutron 
scattering peak was observed in the superconducting 
state of YBa2Cu30j; at a 2D wave vector Q = (ir/a, n/b) 
(a and b are the lattice spacings of the CuC>2 plane) 
and energy u> equal to 41 meV at x ~ 7 |ffl]. This ef- 
fect was confirmed by several groups for x — 7 and 
6.5 < x < 7 The peak has the following remarkable 

features: (a) It appears only below the superconducting 
transition temperature T c f|; (b) It is localized in both 
energy and wave vector; (c) It has sinusoidal dependence 
on q z , the wave vector perpendicular to the CuC>2 planes, 
which implies perfect antiferromagnetic correlation of the 
two planes in a CU2O4 bilayer. 

A number of theories, suggested to explain this effect, 
presume that the peak occurs due to spin-flip electron 
excitations across the superconducting gap ]l|,p|,^[H15[ . It 
was emphasized in Ref. || that, since magnetic scattering 
is odd with respect to the time reversal, the BCS coher- 
ence factor in the neutron scattering amplitude vanishes 
unless Ak has opposite signs for the electron wave vectors 
k and k + q connected by the 2D neutron wave vector 
transfer q (see the inset to Fig. |j]): 



the framework of weakly |7|,|17|,|18[ or strongly correlated 
electrons [[H)] |2l). In this state, the order parameter has 
s-wave symmetry and opposite signs in the bonding and 
antibonding bands. We show below that this state (un- 
like the d-wave state) provides the best explanation of 
the neutron scattering experiments, particularly, of fea- 
ture (c). 

The neutron scattering cross-section is proportional 
to the imaginary part of the electron spin susceptibil- 
ity x(<1j <7zj w )- I n a bilayer system, it is given by the 
following expression p,E0[: 



X (q, q z ,u) = x {+) cos 2 M/2) + X ( - ] sm 2 (q z d/2), (2) 

where d is the distance between the Cu02 planes in the 
bilayer, and 
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Through this condition, neutron scattering can probe the 
symmetry of the superconducting state. For q = Q, 
condition ([!]) is not satisfied for a simple s-wave state, but 
is satisfied for the d x 2_ y 2 state. So, it was suggested in 
Refs. HH-fD that the peak in question is a manifestation 
of the d x 2_ y 2 pairing in YBa2Cu307. 

However, since most of these theories, as well as Ref. 
Jl6| |, dealt with only one CuC>2 layer, they were unable 
to consider important feature (c). Because YBa2Cu3C>7 
consists of CU2O4 bilayers, there should be two electron 
bands formed by bonding and antibonding states. An im- 
portant issue is the relative sign of the superconducting 
order parameter in these two bands. Along with a regular 
d-wave state which has the same sign in the two bands, 
another state, which we will call s ± , was discussed in 



Here, b and a refer to the bonding and antibonding elec- 
tron bands of the bilayer. The susceptibilities x ac- 
count for the transitions between the respective bands. 

In experiments jjjH]], only the second term in Eq. (j|) 
was observed (feature (c)). According to Eq. (||), the co- 
efficient x^ + \ which appears in the first term of Eq. (g), 
involves transitions within the same band. Following the 
coherence factor arguments, we conclude that condition 
(|l|) is not satisfied within the same band, that is, the gap 
has s-symmetry. On the other hand, the coefficient 
of the second term of Eq. (g), which involves transitions 
between the different bands (Q), is not suppressed. Thus, 
the order parameters of the different bands, A( a ) and 
A^ b \ have the opposite signs. This means that the state 
is the s^ state described above. In this picture, the neu- 
tron scattering peak is due to the excitation of electrons, 
say, from the bonding band below the superconducting 
gap to the antibonding band above the superconducting 
gap, and the gaps have opposite signs. 

In order to illustrate the above discussion, we per- 
formed explicit calculations for a well-known realistic 
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tight binding model of YBa2Cu3C>7, which has the fol- 
lowing electron dispersion law: £k = — 2t(coa(k x a) + 
cos(k y aj) — At' cos(k x a) cos(k y a) — fx, where t = 250 meV, 
and t'/t = —0.45 p^j . We have chosen the Fermi energy 
/i = —440 meV, so that the van Hove singularity lies at 
an energy £ v h = 10 meV below the Fermi level p2| , p3[ . 
The corresponding Fermi surface is shown in the inset to 
Fig. [j]. To simplify the calculations, we set the hopping 
amplitude between the layers, t±, equal to zero. This as- 
sumption does not qualitatively change our conclusions. 

In BCS theory, at zero temperature, the susceptibilities 
are given by the following formula pij : 
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where E^ = v^k+^k i s tnc quasiparticle dispersion 
law in the superconducting state, T is the damping con- 
stant, and the indices i and j label the bonding and an- 
tibonding electron bands. Because of the simplifying as- 
sumption t± = 0, the dispersion laws of the bonding and 
antibonding bands are identical and do not need to be 
distinguished in Eq. (||) . Only the signs of the order pa- 
rameters may depend on the index i. Note that only 
when ([j]) is satisfied, does the coherence factor (the first 
line of Eq. (||)) not vanish at small energies £ <C A. 

We have calculated the spin susceptibility Xo directly 
from Eqs. (g) and (fi) using T = 1 meV and the gap 
computed in Ref. E7] for T c — 90 K. The gap attains its 
maximal value Ao = 17.5 meV near the points X and Y 
in the Brillouin zone (see the inset to Fig. [TJ) . Imaginary 

Xq (Q,oj) and real Xo^Qi^) parts of the suscepti- 
bility are shown in Fig. [l] for the superconducting and 
normal states for two values of the Fermi energy p. We 
observe that in the superconducting state, for the en- 
ergies ijj lower than the absorption threshold E t h ~ 35 
meV ("spin gap"), Xo( w ) i s zer o. Furthermore, for the 
realistic /i = —440 meV both real and imaginary parts 
(solid curves) have sharp peaks, at 35 meV and 38 meV, 
respectively. In order to clarify the origin of these fea- 
tures, we repeated the calculations (dashed curves) with 
an unrealistic choice of the Fermi energy /i = —370 meV, 
which moves the van Hove singularity much deeper below 
the Fermi level: £ v h = 80 meV. In this case, the peak in 
Xo(w) stays at the same energy E t h, and xd't^) develops 
a step at the same energy. On the other hand, the peak 
in XoO^O shifts to much higher energy of about 100 meV, 
which is close to A + £ v h- xbfa) develops a negative 
step at the same energy. Similar results were obtained in 
Ref. [E3 for the d x 2_ y 2 state. 

To gain a qualitative understanding of this behavior, 
we set T — > and set the coherence factor in Eq. (|s|) to 
1. In this approximation, Xo(cIj w ) 1S proportional to the 
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FIG. 1. The imaginary Xo' ^Qi^) an d real Xo' ^Qi^) 
parts of the electron spin susceptibility in the superconduct- 
ing s (solid and dashed lines) and normal (dots) states at 
T = 0. Dashed lines and dots correspond to the Fermi energy 
fi = -370 meV, solid lines to fi = —440 meV. Inset: The 
Fermi surface for [i — —440 meV. 



joint density of states A(q,uj) = ^2 k 5(u> — i?k+q — ^k)- 
The two- particle energy ^(k, q) = i?k+ q + E^, consid- 
ered as a function of the 2D wave vector k for a fixed 
q, has a minimum and several saddle points. The min- 
imum defines the threshold energy E t h- In a 2D case, 
the joint density of states has a step at the threshold 
and logarithmic divergences at the saddle point energies. 
Correspondingly, the real part, Xo( w )' which is related 
with Xo( w ) by the Kramers-Kronig relations, has a loga- 
rithmic singularity at the threshold. Exactly that behav- 
ior is observed in Fig. [|. The minimum of i?2(k, Q) is 
achieved at a vector k such that both k and k+Q belong 
to the Fermi surface (see the inset to Fig. [j]), where both 
i?k and -Ek+Q attain their minimal values approximately 
equal to A . Thus, E th « 2A = 38 meV. 

The logarithmic peak in Xo( a -0 occurs because of tran- 
sitions between the occupied states located near X and 
Y points and empty quasiparticle states above the su- 
perconducting gap. The points X and Y are the saddle 
points of the normal-state dispersion law £k- They pro- 
duce the van Hove singularity in the single-particle den- 
sity of states in the normal state at the energy £ v h- The 
logarithmic divergence in the joint density of states in 
the superconducting state is located, then, at the energy 
E* k A + + £vH- The value of £ v h depends on 

the Fermi energy pi. For the two choices of /i in Fig. [I], 
the values of £ v h are equal, respectively, to 10 and 80 
meV, which gives values of E* w 38 meV and E* ps 100 
meV, respectively, in agreement with the positions of the 
peaks in Xo( w ) m Fig. [j]. Note that both LDA calcula- 
tions and photoemission experiments [^3| place the 
van Hove singularity very close to the Fermi level, which 
corresponds to the first choice. 

Comparison of Xo(Q,^>) m the normal state at zero 
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temperature (dotted curve in Fig. [T]) and in the super- 
conducting state (dashed line) shows that the logarithmic 
peak due to the van Hove singularity is absent in the nor- 
mal state. Instead, we see a cusp at an energy close to 
£ v h = 80 meV, which marks the threshold where transi- 
tions from the saddle points to the Fermi level become 
allowed. The reason that, instead of a cusp in the normal 
state, a divergence develops in the superconducting state 
is that in the latter case, roughly speaking, the transi- 
tions take place between the van Hove singularity and 
the coherence peak in the quasiparticle density of states. 
Since there is no threshold of absorption in the normal 
state, Xo( w ) nas n0 divergence at a finite us in this state 
(dotted line in Fig. |l|). 

Summarizing, we conclude that the peaks in Fig. [I] ap- 
pear to manifest two different physical mechanisms. The 
peak in Xo(Q> w ) ^ s °-ue to the transitions between the 
states just below the gap to the states right above it at 
the energy E t h ~ 2Ao which depends only on the super- 
conducting gap. On the other hand, the peak in x'o (Q, w ) 
is due to the transitions from an occupied saddle point 
to the empty states just above the superconducting gap. 
The energy of this peak depends on the van Hove en- 
ergy of the normal state £ v h- For a realistic Fermi en- 
ergy fx = 440 meV, £ v h = 10 meV is small compared to 
A = 17.5 meV, therefore the peaks in x'o and Xo ap- 
proach each other, which is shown below to have a major 
effect on neutron scattering. None of these peaks appear 
in the normal state. 

Curves in Fig. [I] are reminiscent of the experimental 
curves where x"(Q,w) was found to exhibit a spin 
gap and a peak at different energies in underdoped sam- 
ples with x<7. It is tempting jl4| to identify the log- 
arithmic peak in XoiQ^) with the peak observed in 
the experiment. However, the plot of Xofa u ) reveals a 
structure quite different from that found experimentally. 
While there is a local maximum at q = Q and u> « 2Aq 
in Fig. ^, the corresponding peaks, although weaker, ex- 
ist also at other values of q, in contradiction with feature 
(b). Even worse, when q^ 0, Xo(q,o;) increases to much 
higher values than at q = Q. We conclude that x'o (Qi w ) 
does not provide a satisfactory explanation of the exper- 
imentally observed neutron scattering. 

Thus far, we have neglected the interaction between 
quasiparticles in the superconducting state. If we take 
into account a weak antiferromagnetic interaction be- 
tween Cu spins in the plane, J(q), and between the 

planes, J±, the bare spin susceptibilities Xo \ dis- 
cussed above, should be renormalizcd in an RPA manner 
0[||OU2|: 

X (±) (q,^) = *£±0(q,a;)/[l + j(±)(q) x ( ± )(q 1 a,)/2], (6) 

J (±) (q) = J(q)±J±. (7) 

Since the real part, x'o Hht^j diverges as us — + 
E t h, as discussed above, the renormalized susceptibility 




FIG. 2. Bare susceptibility x'o ilx , q y , ui) , plotted as a 
function of q — q x = q y and us for the s order parameter 
with (U = —440 meV. 

(q, us) (^J) has a pole at us close to E t h with singulari- 
ties in both real and imaginary parts of X ■ Physically, 
the pole in ^"-'(q, k>) describes a triplet electron-hole 
collective mode (an "antiparamagnon" or a "triplet exci- 
ton") with the energy slightly below E t h- In this picture, 
the peak in the neutron scattering rate occurs due to 
inelastic excitation of this collective mode by neutrons. 

Since the singularity in x'o (Qj w ) exists, generally 
speaking, for any q, the exciton state (and, thus, a peak 
in neutron scattering) should exist for any q, in con- 
tradiction with feature (b). However, if we take into 
account a finite lifetime of quasiparticles T in Eq. (|^), 
the excitons also acquire a finite lifetime. For the an- 
tiferromagnetic interaction between the nearest neigh- 
bors, x acquires a peak at q = Q, where 
J(q) = Jii [cos^a) + cos(q y b)] reaches its maximal neg- 
ative value. Thus, the position of the neutron peak in 
q-space is set by J(q) and in us by Xo (Q)'*')- This 
statement is illustrated in Fig. [|, where we show x"(q> w ) 
calculated according to Eq. (g) with — 150 meV and 
r = 1 meV for the state. In agreement with exper- 
iment, a single peak in X w )) localized both in q 
and us, is observed, which is now more than twice higher 
than that at q — > 0. The magnitude of the peak depends 
on the chosen value of J and T, but the qualitatively the 
picture remains the same for a reasonable range of these 
parameters. Similar results were obtained in Ref. jTq ] 
for the d x 2_ y 2 state. It is worth noting that the peak in 

X is due to the peak in the real part of Xo ■ When 
the van Hove singularity is close to the Fermi level, the 
peak in x'o ' (Qj u ) gets stronger and enhances the peak 
in the renormalized susceptibility X (Qj^)- 

The crucial difference between the and d x 2_ y 2 states 
is due to the following: In the state, the divergence oc- 
curs only in x'o i but not in Xo + ' where it is suppressed 
by the coherence factor (pt) . This suppression, unlike the 
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J =150 meV 




FIG. 3. Renormalized susceptibility \' ~ (ix, q y ,^>), plot- 
ted as a function of q — q x = q y and uj for the same s ± state 
as in Fig. H. 



case of logarithmic divergence in Xo 1 i s exact, because 
the divergence in the real part of Xo + ^ would occur pre- 
cisely at the threshold energy where the coherence factor 
vanishes. As a consequence, among the renormalized sus- 
ceptibilities, only ^ diverges, but not X ■ Taking in 
account Eq. (^), this naturally explains feature (c) of the 
experiment. 

In the case of the d-wave, the coherence factor allows 
divergences in both Xo ' an d Xo + "* > which, in turn, pro- 
duce divergences in both x an d X ■ The only way 
to reconcile this with the experiment is to assume that 
is small: | J^ \ -C or even has a wrong (pos- 

itive) sign. According to Eq. (^), this would require that 
the antiferromagnetic interaction between the layers be 
stronger than within a layer: J± > 2Jii, which is unreal- 
istic §§. 

A number of theoretical papers J7j-pd[| consider the spe- 
cial case when t' = in the dispersion law. This as- 
sumption results in nesting at the Fermi energy or at 
another energy ("dynamic nesting"), which produces a 
peak in x '(Q,w). However, this case is not relevant for 
YBa2Cu3C>7 where Q is not a nesting vector. 

In conclusion, we considered a scenario where the peak 
at 41 meV observed in neutron scattering experiments 
[Q-Q occurs due to spin-flip interband electron excita- 
tions across the superconducting gap which are enhanced 
by the antiferromagnetic interaction between Cu spins. 
We found that the experiment can be explained most 
naturally if the superconducting order parameter is of 
the type, that is, the order parameter which has the 
s-wave symmetry and the opposite signs in the bonding 
and antibonding bands. This state easily explains the 
observed dependence of the scattering intensity on the 
momentum perpendicular to the Cu02 planes. On the 
other hand, the d x 2_ y i case can be reconciled with the 
observed dependence only if the antiferromagnetic inter- 
action between the Cu02 planes is stronger than the in- 



teraction within the plane, which seems to be unrealistic. 
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